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I. INTRODUCTION 



For a long time, much attention has been concentrated on the time-dependent harmònic 
oscillator both in classical and quantum mechanics. The quantization of a time-dependent 
harmònic oscillator is very important when we treat free quantum fields in curved space- 
time M, quantum statistical properties of radiation [FJ, coherent states and squeezed states 
of light §, etc. 

One of the most powerful methods to find a quantum state of a time-dependent oscil- 
lator is the generalized invariant method introduced by Lewis M and Lewis and Riesenfeld 
(LR) ||. They have derived a simple relation between eigenstates of LR invariant and Solu- 
tions of the corresponding Schròdinger equation. The Schrodinger wave functions have been 
studied in Refs. ||[7j and their results are generalized in Ref. [§]. 

The time-dependent harmònic oscillator with damping and a perturbative force was 
considered by Khandekar and Lawande who found the exact wave function in terms of an 
auxiliary function p and an external force f(t) by using the LR invariant |J. The exact 
evolution operator for the harmònic oscillator subject to an arbitrary force was obtained 
by Lo ||10|| . He obtained the quantum states by solving the differential equation satisfied 
by the evolution operator. For a generalized harmònic oscillator (GHO) with a constant 
force, the invariant formulation was developed and the geomètric phase was calculated in 
Ref. | Ï2|| . Recently, the authors have calculated the generalized invariant and its eigenstates 



for a time-dependent forced harmònic oscillator in the Heisenberg picture [|TTJ. In addition, 
they compared their quantum states with them of Lo. 

In this paper, we extend our previous work [[ÏÏ]] to a general quadratic oscillator. By 
adding time-dependent terms proportional to pq + qp and p to the Hamiltonian of Ref. (TT| 



we have the most general type of a time-dependent quadratic oscillator. Hereafter, we call 
this type of oscillator a generalized driven oscillator (GDO) to discriminate it from a GHO 
which has no linear driving terms. This work is also an extension of Ref. in the points 
that a driven force is generalized from the constant one to a time-dependent one and the 



2 



p term is not considered in Ref. [0. For a GDO, we find the LR invariant, its eigenstates 
and the corresponding Schrodinger wave functions. Analyzing the wave functions, we find 
the CISs and calculate the corresponding non-adiabatic Berry phases. 

In Sec. y we find a generalized invariant for the GDO. First, we find a generalized 
invariant for a GHO, where we represent the invariant in terms of the classical solutions. For 
this invariant, we construct the creation and annihilation operators. Using these operators, 



we find a generalized invariant for the GDO employing the method introduced in Ref. [JTT 



We also construct the Fock space using the eigenstates of the invariant. In Sec. |TJ we obtain 
the exact wave function which satisfies the time-dependent Schrodinger equation. In Sec. IV 
we study the condition for the existence of CISs and calculate non-adiabatic Berry phases. 
In Sec. [V] some wave functions are presented as examples. The last section is devoted to the 
discussion. 

II. GENERALIZED INVARIANT AND ITS EIGENSTATES 

Consider the Hamiltonian of a GDO: 

H T (t) = H(t) - F{t)q(t) - G(t)p(t), (2.1) 

where 

H(t) = ^/(t) + ^ W)q{t) + q(t)p(t)\ + l -M{tW{t)q\t) (2.2) 

and M(t),u>(t) 2 , and Y(t) are time dependent parameters that satisfy Lü 2 (t) — Y 2 (t) > 0. 
The F(t) and G(t) terms describe linear driving mechanism. 

Since the set of Hermitian operators, {^p 2 (t), |ç 2 (í), \[p{t)q(t) + q(t)p(t)],p(t), q(t), 1}, 
forms the Lie àlgebra, we try to find the following form of the LR invariant: 

I T (t) = I(t) + gi (t)q(t) + g 2 (t)p(t) + g 3 (t) (2.3) 

where 

n*) =9-{t)—=- +go{t) 7: + g + {t)— — . (2.4) 
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A. Generalized Harmònic oscillator with no linear driving terms 



Before proceeding with the Hamiltonian ( |2.1|) , let us consider the Hamiltonian of the 



GHO (|2.2|) . The invariant for this Hamiltonian will be the form of (2^), and from the 
invariant equation (in % = 1 units) 

p(t) = ^I(t)-i[I(t),H(t)] = Q, (2.5) 

we obtain a coupled linear system of the first-order differential equations 

2 

9- = 2 9- Y ~ jj9o, 

g = Mu?g_ - —g + . (2.6) 
g+ = -2g + Y + 2Mu 2 g , 

where a dot denotes a time-derivative. Introducing 

9- = a g-, g = 9o, 9+ = — , M(t) = — — , 2.7 



where a(t) = e~ J Yd \ it follows from (|2~^) that 
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9- 




fio 


= Muj 2 g- - 


h 


= 2Muj 2 g . 



1 . 



(2.8) 



Here we note that 9±,o satisfy the same differential equations as Eq. (6) in Ref. |Ï3| . Thus ^±,0 
can be represented in terms the solutions of the differential equation (d / dt)[M f]+ Mu 2 f = 0. 
Further, we have the relations fi = afijfy = a/2, where fi^if) are the two independent 
solutions of the equation of motion for H(t): 

Mft + Mfi + (Meu 2 - MY 2 - MY - MY)f { = 0. (2.9) 

Therefore the most general solutions of Eq. (|2.8| ) are given by 
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g.(t) = c t fl{t) + c*h{t)h{t) + c 3 / 2 2 (t), 



,/„(/) = ( dA(t)/i(t) + | [/i(t)/ a (t) + /l(t)A(t)] + C 3 f 2 (t)j 2 (t) ) - (/),/_(/). 



<7 + (t) = M 2 (í) (c^ (í) + c 2 h(t)f 2 (t) + c 3 / 2 (í) ) + M(Í)F 2 (Í)^(Í) + 2Y(t)g (t), (2.10) 

where ci,C2 and C3 are arbitrary constants. 
Now the invariant (12.41) can be written as 



I(t)= Ul (òt(í)ò(í) + i 



introducing 



ò(í) 
6t(t) 



2g_(t) y2üj ig -(t) 



' ffo(í)U(í)+M/^p(í), 



2u, 



s y2g-(t) V 2üJ ig „{ty^~' )^~' 1200! 
with o?j = g+(f)g-(t) — g%(t). From the equation of motion for 6(í): 

>^K Í )-,K Í ),,(^-^K Í ), 

we have 



(2.11) 



(2.12) 



(2.13) 



6(í) = ò e - ie(í) , 



(2.14) 



where 



e(t) 



dí 7 - 



Ito M(t')g-(t>) 

Then it is clear from ( 2.11|) and (|2.14 ) that I(t) is an invariant for H(t). 



(2.15) 



B. General quadratic Hamiltonian with linear driving 



Now we find the invariant ( |2.3j ) for the total Hamiltonian ( |2.1| ). Inserting the ( |2.1D and 
3|) into 



|/t(í) = |Jt(í)-í[/t(í),ÍÏt(í)] = 0, 



(2.16) 



we get the linear coupled differential equations ( |2.6|) and additionally we have for «71,2,3 (£): 



<7i = -Y 9l + Mu 2 g 2 + {Gg + - Fg ), 
92 = ~T79i + Yg 2 - (Fg_ - Gg ), 
9a = Ggi - Fg 2 . 



(2.17) 



However these equations are too difEcult to get the solutions, and we get around this difficulty 
using the method introduced in Ref. [II] . 

For the presence of the linear driving terms we search for the invariant of the form 



I T {t)=uj I (B\t)B{t) + - 



(2.18) 



where 



B(t) = b{t)+P(t). 



(2.19) 



It follows from 



j t b(t) = ^b(t)-i[b(t),H T (t)} 



—1 



M(t)g_(t) 



b(t)+W(t), 



(2.20) 



with 



W(t) 



+iJ^^9o(t))G(t)+J 9 -^F(t) 



2g_(t) \2uj ig _(ty 



2u, 



(2.21) 



that 



ÉL B (t) = -i—^L—B(t), 
dt V ; M(t)g-(t) V ; 

provi ded that f3{t) satisfies the differential equation 



(2.22) 



(2.23) 



The solutions of Q2.22 ) and ( |2.23| ) are easily found: 



B(t) = e- i0{t) B(t o ), B\t) = e ie ®BÏ(t ), 



(2.24) 



and 



(3(t) = e- i0 W (fa - f W{t')e ie ^dt^ 



(2.25) 



where (3q is an integration constant. 

Now it is trivial to find g t {t) (i = 1,2,3), by substituting (|2~T2j ) into ( gTg ) using (£39]) 

we get 



J T (Í) = I(t) + vfaF^rq® + UjJ^ í p(í) + 



+ ^(^ + /3 7 2 ), 



(2.26) 



where /3r (Pi) is the real (imaginary) part of /?. By comparing 
easily read that 



gi(t) = lüiPrJ + ujiP 



g 2 (t) = lútfi 



9-(t) 

l ig-(t) 



l ig-(t) g Q (t) 

íúi g_(t)' 



3[) and Q2.26j ), we can 



(2.27) 



g 3 (t)=u I (Pj i + (3j). 



(2.28) 
(2.29) 



Thus we have found the generalized invariant - ( [2.3|) or (|2.18|) - for the GDO ( [2.1|) . 

As in the usual case, we construct the Fock space as eigenstates of Ir(t), which are given 

by 

£t n (í) 



\n,t) 



\Q,t) 



(2.30) 



where 



B(t) |0,t) = 0. 



The eigenvalue equation of Ir(t) can be written as 

I T (t) \n, t) = uji {n + ^\ \n, t) , n = 0, 1, 2, 3, 



(2.31) 



(2.32) 
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III. WAVE FUNCTION OF THE SCHRODINGER EQUATION 



Following LR H the wave function of the Schròdinger equation 



i^\Mt))s = H T (t)\Mt))s (3-1) 



can be written by 

\Mt)) s = e ta " {t) \n,t), (3.2) 
where the phase functions a n (t) are found from the equation 

ja n (t) = (n, t\ ^ - \n, t) . (3.3) 

A. The coordinate representation of the eigenstates of the invariant 

Let us define new coordinate Q and momentum P by taking the following time- dependent 
successive unitary transformations: 

Q(t) = Ul(t)q(t)U 1 (t)=q(t), (3.4) 

P(t) = Ul(t)p(t)Ui(t) = Pit) + ^L(t), (3.5) 

Q(t) = U\{t)Q{t)U 2 {t) = Q(t) + A q (jt), (3.6) 

P(t) = U%(t)P(t)U 2 (t) = P(t) + A p (í), (3.7) 

Q(t) = U}(t)Q(t)U 3 (t) = - T L=Q(t), (3.8) 

P(t) = Ul(t)P(t)U 3 (t) = y/Ï4£jP(t), (3.9) 



where 



U x {t) = exp(i^A (3.10) 
U 2 (t) = exp (-iA q (t)P) exp (iA p (t)Q), (3.11) 
C/ 3 (í) = exp Q(PQ + QP) ln^(í)) , (3.12) 



with 



Then the LR invariant can be rewritten as a simple form 



The eigenfunction of Ir{t) in Q coordinate is well known as 



Qi t\n 



íúi 



1/4 



ui Q 

e '~H n [y/ÜJïQ) 



V2™rï! V 7r 

where iï n is a Hermite polynomial and which is normalized to be (n, t\n',t) = 5 nn r. 
Now we represent the wave function in terms of original coordinate variable 

(q\n,t) = <p n (q,t) = (q,t\n) . 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



It follow from Eqs. (PHQ , ( pi]) , and ( ^T2|) that 



e l °- 9 {Q,t\ 



e -if-q*-iQA p 



(9- 



-1/4 



Q,t 



\9-J 

and therefore the coordinate representation of the invariant eigenstate is found to be 



(3.17) 



<Pn{q,t) = (q\n,t) 



(q + A g ) . (3.18) 



We can show that Eq. (|3.18| ) satisfies the equation lT(t)(p n (q, t) = Lüi(n+^)(p n (q, t) by direct 
differentiations. It should be noted that the coordinate representation of the eigenstate is 
not unique since the time-dependent phase is arbitrary. Consequently, we should find a n (t) 
in (|3.2|) by solving Eq. (|3.3|) to find the wave function satisfying the Schròdinger equation. 
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B. Schròdinger Wave Function 



As in Ref. H, we find the phase a n (t) in Eq. (|ü.2p using the following relations 



[n, t\ ^- n, t) = {n- 1, í — /i - 1, í) + -= (n, t\ — — \n, t) 
ot ot \/n ot 



«9í W M(t)g.(ty 



B\t)+i\B\t),H T (t) 



and 



F T (Í) = # Q (í) - B(t) (h-(t)/3(t) + -h (t)(3\t) + iW* 



B\t) ( /i+(í)/3+(í) + -h (t)(3(t) - iW(t) 



+ (iW*(t)P(t) - iWft(t)) + \ (h^(t)(3 2 (t) + M(í)^(í) + M*)/^ 



where 



and 



/i (í) 



i 



M(í)(?_(í)^ 
1 



<? 2 (í) + M\t)u\t)gí{t) + W/ 2 - 2M(t)Y(t)g (t)g_(t) 
g 2 Q (t) + M 2 (t)iü 2 (t)g 2 _(t) - oü 2 j T 2ig (t)uj 



2M(t)g_{t)uj I 
-2M(t)Y(t)go(t)g-(t)±i2M(t)Y(t)g-(t)Lü I ] . 



After a little àlgebra, we obtain the phase 

oí n vt) — —{n-\ — ) / dt — — — — 
w v 2 J J M(t')g_(t') 



+ fdf [o£(ío 



Ui 



2üü! 



M(t')g4ï) V^-(í') 
Therefore, we find the exact solution of the Schròdinger equation: 



X 



V2^.\ng-J Vv 9- 



This wave function becomes the one of Ref. [[y| when F(t) = = G(t) and Y(t) = í 
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IV. EXISTENCE OF CYCLIC INITIAL STATES AND NONADIABATIC BERRY 

PHASES 



In this section we examine carefully the wave function (3.25) to find CISs and the cor 



responding non-adiabatic Berry phases. Suppose that all parameters M(t), Y(t), u(t), F(t) 
and G(t) are periòdic function of t so that the Hamiltonian is cyclic with a period r. The 
definition of CIS gives ip(t Q + r) = e* x ?/>(ío), and the Berry phase T is given by T = x — 8, 
where 5 = - % (ip(t)\H(t)\ip(t))dt. 

Inspecting closely the form of the wave function ( |3.25|) , we find that the necessary and 



sufficient condition for the existence of CISs is the existence of periòdic g~(t) and /3(t). The 
periòdic property of g~(t) is related to the solution of Eq. ( |2.9| ). The conditions for the 



periòdic solutions of classical equation of motion were well studied in Ref. ||13|| . Therefore 
we investigate only the condition for the periòdic (3{t). Hereafter we assume that g-(t) is a 
r-periodic function. 

To study the periodicity of fl(t), from (|2.25|) we write 

rt+T 



P(t + t) = e - t0{t+T) (A) - / T W(t')e ie ^dÍ 

'(3(t) - aoe-^W] , (4.1) 



e -i6o 



where we used Eq. ( ^.25[ ) again. Here 9q and o"o are defined as 



dt' , , : (4.2) 
and 

' +T W(t')e ie{t,) dt'. (4.3) 



0o 



As seen from the periodicity of the integrand of ( |4.2j ), 6 is constant. Then, if e~ l0 ° = 1, the 
periodicity of F(t) and G(t) guarantees the periodicity of the integrand of (^3|) [see ( |2 . 2 1| ) ] , 
and hence <To is also constant. Then it is clear from ( |4.1| ) that the conditions for periòdic 
(3{t) are (i) 9q = 2im(n = integers) (ií) o"o — 0. When these conditions are fulfilled, all the 
eigenstates of the generalized invariant are CISs with 

11 



Xn = a n (t + r) - a n (t) = a n (r) - a n (0) 



(4.4) 



and the nonadiabatic Berry phase is obtained by 



r„ = xn+ / (Mt)\H(tM n (t))dt 

Jo 



Using ( |3.21| ), we get 



r„ r = n + ~ 



1 



dt 



+ / dt 
Jo 



2/ jo 
~ho(t) 



9 2 o(t) 



Y(t)g (t) 



M(t)g-(t)üJi uo 



where £r is the real part of £ (£) with 



$ = 2i r(t) 



(4.5) 



(4.6) 



(4.7) 



C= M*) + 



/? 2 (t). 



(4.8) 



M(t)^_(í) 

Here we have taken partial integration using Eq. fl2.6p to get the first term of (|4.6|) as in 
Ref. When -F(t), G(í) and K(í) vanish, this Berry phase reduces to the one of Ref. JT3 



Note that the minus sign of Eq. (4.8) [Eq. (5.2)] of Ref. |y§ is an error. This error can be 
easily checked from another expression for the Berry phase Eq. (4.7) there. 



V. EXAMPLES 



In this section we consider some physically interesting examples and write the explicit 
forms of their wave functions and study on their Berry phases. 

Example A. M(t) = m, u(t) = ujq, F(t) = Fq, and Y = = G: The two independent 
solutions of (|2.9|) are chosen to be fi(t) = e %u)ot and /2(í) = e~ luot . We can set g-(t) — 1/m, 
lüj = üüq by setting c\ = = c 3 , c 2 = 1/m in (|2.10|) . Further, it follows from ( p.23|) that 



at v ImujQ 



(5.1) 



and its general solution is given by (|2.25|) as 
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P = **** Po + n 7^= (5-2) 



It is straightforward to find the wave function from ( |3.25 ), however we confine ourselves to 
the special case, where f3(t) = — l· , ^ mJ0 ■ Then the wave function is given by 



xe 



o H n (y/mwo (q - F /mc^)) , (5.3) 



which shows clearly that the wave function of the oscillator is shifted by the amount of 
F /muj 2 . We can always take an arbitrary period r, then 6q = Uot = 2j ^ k , i.e. r = 
Further, cr = 0. Thus all the eigenfunctions of the invariant of this system are CISs, but all 
the Berry phases vanish as expected from the time-independency of the system. 

Example B. The Caldirola-Kanai oscillator |Tj|ïü]: The Hamiltonian of the Caldirola- 
Kanai oscillator with external force is given by 

"W-i^ + ^pT^ + H· (5 ' 4) 

which describes, at classical level, a oscillator with a time-dependent frequency and a velocity 
dependent damping term. The equation of motion without the external force is given by 

q(t) + 2jq(t) + u 2 q(t) = (5.5) 

whose solutions are 

fi(t)=e^ t+int , f 2 (t)=e-*- iSH , (5.6) 

where Q = ^/u 2 — 7 2 . By setting c\ = = c 3 and c 2 = 1/m we get 

g-{t) = -e-^\ flb(t)= 7 , 9 + it)=m^e 2 ^. (5.7) 
m 

Then, the wave function is given by 



13 



x exp 



( v^e* ( q + J—e^Qn ] ] , (5.8) 



where the function (3(t) is the solution of the following equation 



and its solution is given by (|2.25|) . Unless 7 = 0, the function g~(t) is not periòdic and this 
system has no CIS. 

Example C. Special case of Example B with 7 = 0, /(£) = — F sin(o; e í), with (u e 7^ 0). 
In this case we have 

9-{t) = -, 9o{t) = 0, g+(t) = mu\ (5.10) 
m 



and 



d „, x . ~, x .F sin(o; e í) 
The solution /3(í) is given by 

Tp / p ÍLü e t p — Ííüt p ~ÍLü e t p — ÍCJÍ \ 

P(t) = fre** + i-^= % e - —, 6 — , (5.12) 

for tu 7^ u e . For /3(í) to be periòdic, it is sufficient that uj/uj e is a rational number, i.e. 
uj/uj e = r/r e , where r and r e are positive integers not commensurate. Then the system has 
the CISs with a period r = 2nr/uj = 27rr e /u e . One note that the period of this system is 
r e = 2ir/u e , while the CISs are r e r e -periodic. Using h = 2uj and /i_ = 0, from (|4.6|), the 
Berry phase is given by adding the following three terms: 







T h (t)\P(t)\ 2 dt = 2nr 



2 F ((3 - (3* )u e F 2 ( u' + ul 1 



tV2^(uo 2 - ui) 2muj - ulf 2{uj 2 - cu 2 ), 

(5.13) 
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and 



F 2 1 

Ç R dt = -litr-^ -, 

o Zmuj uj z — ují 



r T F 2 1 
/ ( R dt = 27rr-^ -, 



(5.14) 



(5.15) 



as 



Y n = 2-kt 



\f%\ + 



F (A)-/3 >e 



+ 



F 2 

^0 



+ w 2 



(5.16) 



iy/2muj{uj 2 - cj 2 ) 2muj \(üj 2 - lú 2 ) 2 u 2 - u 2 _ 
It should be noted that, for this example of a forced harmònic oscillator, the Berry phase 
of the eigenstate \n,t) is independent of quantum number n [in Eq. ( |4.6|) the first integral 
term which includes n, vanishes]. 



Example D. The damped pulsating oscillator 10 



M(t) = m e 2 ^ t+l·lsinut \ u 2 (t) =íl 2 + , 1 —jM(t). (5.17) 
The motions of quantum operators are analyzed in Ref . jïï] . Now we find the wave function 



using the classical solutions and the auxiliary functions used in Ref. |11 



í .1 d 



e ■§ Tí + § M sin i/í 



\j2 n n\ V 7T / 
(t)q 2 - A p q - Knoüe 2 ^ 2 ^ 1 (q + A q )*j 

xH n ( yjm ne^ + ' taavt (q + A q )) , (5.18) 



x exp —i——m 
F V Adt 



where 



A, 



-■yt—fi sin vt i 



fa, A p = ^moSle^^fa. 



(5.19) 



Here j3{t) satisfies 



(5.20) 



and the solution is given by ( |2.25| ) . As in the case of Example B this system cannot possess 
a CIS because of the damping effect. 
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Before closing this section we find the variances in p and q for a quantum state \n,t). 
From (gTg ) and ( |2TT9|) we find 



( n |[Aç(t)] 2 |n) = (2n + l)^, (5.21) 



(n|[Ap(í)] 2 |n) = (2n + l) 



(5.22) 



2»-(í) 

We note that as in Ref. |TT[ the linear driving terms don't change the width of a wave packet 



but change only the position of the wave packet. 

VI. DISCUSSION 

For the most general quadratic harmònic oscillator, we obtained the generalized invariant 
It and its eigenstates. For these eigenstates, we have obtained the Schròdinger wave function 



in a closed form, which can be applied to finding the wave functions of our previous work |] I 
The particular point is that the quantum wave function can be written as the solutions of 
classical equation of motion with no linear driving terms. The linear driving terms [F(t) 
and G(t)} merely shift the wave packets of the eigenstates of the invariant with no linear 
driving terms by the amount of A q - see ( |3.25| ). We have also investigated on the existence 
of a CIS, and have found the conditions for the existence of CIS. For the CISs, the explicit 
form of Berry phases were obtained. 

The Hamiltonian ü/y (i) studied in this paper, by writing 

can be written as 

H T (t) = Si(í)a 2 + Si(t)a f2 + s 2 {t)[a ] a + aa ] ] + s^{t)a) + s^(t)a, (6.2) 

where 

si = ~Y, s 2 = y s 3 = ~/=j {~ F + imuG) . (6.3) 
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This Hamiltonian is studied in Ref. [|T^] to describe the photon process driven by a classical 
time-dependent source. The time-evolution operator and geomètric phase of this system 
appear in Ref. JÏ9] . Therein the physical quantities are expressed in terms of some auxiliary 
functions. While in our formulation they are expressed in terms of classical solutions of the 
corresponding GHO. Therefore we expect that our formulation is useful in describing the 
quantum optical properties of the model. 
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